% Ejemplo de invocaci??n
% rk4(0,48,10,0.01,'f1');
% como ultimo parametro recibe la funci??n del modelo
% estan disponibles f1 y f2 (el que tiene u)

function energy=Puntod2()
% Parametros de Runge kutta-4
a = 0;
b = 48;
ya = 10;
h = 0.01;
func = 'f2';

% modelo de simulaci??n
kArr=[0.25,0.5,1,2,5,10,15,20];
global ku;
global k;
global TOUT_M;
global TD;

TD=20;
TOUT_M=10;
ku=0.7;
colors = {'1', '2', '3', '4', '5', '1*', '2*', '3*'};
labels = {'k=1/4','k=1/2','k=1','k=2','k=5','k=10','k=15','k=20'};

% Curvas de la temperatura T(t) para distintos valores de k
for i=1:length(kArr)
	k=kArr(i);
	RK = rk4gral(a,b,ya,h,func);
	energy(i,1)=k;
	energy(i,2)=trapz(RK(:,1),RK(:,2));
end

plot(energy(:,1), energy(:,2), '-r*');

%title('Energia para los distintos valores de k (segundo modelo)')
xlabel('k');
ylabel('Energia');

end
